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Abstract. We report on a systematic implementation of su{2) invariance for matrix 
,^ ' product states (MPS) with concrete computations cast in a diagrammatic language. 

As an application we present a variational MPS study of a spin-1/2 quantum chain. For 
efficient computations we make systematic use of the su{2) symmetry at all steps of the 
calculations: (i) the matrix space is set up as a direct sum of irreducible representations, 

jrt . (ii) the local matrices with state-valued entries are set up as superposition of su(2) 

tr. ' singlet operators, (iii) products of operators are evaluated algebraically by making 

use of identities for 3j and 6j symbols. The remaining numerical computations like 
the diagonalization of the associated transfer matrix and the minimization of the 
energy expectation value are done in spaces free of symmetry degeneracies. The energy 
i-rt ■ expectation value is a strict upper bound of the true ground-state energy and yields 

^ ' definite conclusions about the accuracy of DMRG results reported in the literature. 

^ • Furthermore, we present explicit results with accuracy better than 10~^ for nearest- 

and next-nearest neighbour spin correlators and for general dimer-dimer correlators in 
the thermodynamical limit of the spin- 2 Heisenberg chain with frustration. 
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1. Introduction 

We use a systematic implementation of su{2) invariance for matrix-product states 
(MPS), which parallels [H |2], and perform a variational sm(2)-MPS study for the 
frustrated antiferromagnetic spin-| Heisenberg chain. 

The class of MPS is of particular importance for the study of quantum spin 
chains. The MPS appear in different, but closely related ways. Historically, products 
of matrices with entries from a local Hilbert space first appeared as exact ground- 
states (matrix-product ground-states, MPG) for Hamiltonians with special local ground- 
state structure, see [S, H] and developments [5]-[l3]. Second, in the literature on 
integrable systems, vertex-operators ^20j-^22j were introduced for a transfer- matrix 
like construction of ground-states of lattice systems. Third, the density-matrix 
renormalization group (DMRG) [I1]-[I7] algorithms were shown f[E\ [T^ to result in 
states of MPS type. The MPG, vertex-operators and DMRG are all realizations of MPS. 
Consequently, this important class of states attracted strong interest in the quantum 
computation community [23]-|29]. 

In applications of MPS, details may differ strongly. For instance, in MPG and 
DMRG realizations the matrix index space is finite dimensional, whereas in the vertex- 
operator case this space is mostly infinite dimensional. In MPG, the MPS are used 
as an ansatz for an "exact state" for which a "parental Hamiltonian" is to be found 
in subsequent investigations. In DMRG, the MPS appear as variational states for the 
Hamiltonian. In case of the vertex-operator, the full ground-state structure is captured 
at the expense of infinite-dimensional matrices. 

Our investigation is motivated by the need for a computationally most efficient 
scheme for general 5^(2) singlet states of MPS type. The local implementation of Lie 
group invariance has been studied in for instance |33j where 5^(2) invariant MPG and 
parental Hamiltonians were investigated, and particularly in the early work [U |2] where 
also the variational analysis of su{2) invariant MPS for quantum chains and ladders was 
introduced. The paper [1] is well-known for the insight that the finite-system DMRG 
leads to quantum states in MPS form, over which it variationally optimizes. To our 
knowledge the computational scheme of su{2) invariant MPS with arbitrary matrix 
space did not attract the attention the papers [H [2] deserve. Non-Abelian symmetries 
were implemented in DMRG studies in for instance [301 EI]- 

Here we use a variational computation scheme very similar to [HE] based on su{2) 
invariant MPS. However, we reduce the necessary constructions to a minimum without 
taking any reference to DMRG algorithms. Also, we are going to use a diagrammatic 
representation of some of the key objects and relations occurring in the process of 
the evaluation of the norm and the energy expectation value of the MPS. We hope, 
this approach will make the subject as accessible as possible. Our main application 
will be the study of ground-state energies, spin-spin and dimer-dimer correlations of 
the (frustrated) antiferromagnetic spin-1/2 Heisenberg chain with nearest- and next- 
nearest neighbour interactions. For not too strong frustration, this system shows critical 
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behaviour. Still, the variational su{2) invariant MPS give excellent results even for the 
correlation functions. Probably, and in contrast to the expectation expressed in [2], 
the approach is applicable even to odd-legged spin- 1/2 ladders that are not finitely 
correlated. 

Obviously, and following scientific lore, it is important to use all available 
symmetries to find invariant blocks of the transfer matrix as low dimensional as possible 
in order to reduce the computational work involved with the diagonalization procedure. 
Even more important than the economical treatment of the transfer matrix is the 
efficient, nonredundant parameterization of the local building elements, i.e. the matrices 
with entries from the local Hilbert space. It is essential to parameterize these objects 
with as few parameters as possible to reduce the computational time of the minimization 
of the energy expectation value. 

The paper is organized as follows. In Sect. [2] we present a fairly general derivation 
of equilibrium states in the form of MPS and shortly summarize the tensor calculus of 
MPS with emphasis on realizations of symmetries. In Sect. [3] we introduce the su{2) 
invariant local objects based on Wigner's 3j symbols. Here we also introduce the transfer 
matrices and evaluate products of operators by making use of identities involving 3j and 
6j symbols. In Sect. IHwe present explicit results from numerical evaluations of the basic 
formulas derived in the previous section. The results are compared with DMRG data 
of the literature [32] for the frustrated spin-| Heisenberg chain and conclusions about 
the accuracy of the methods are drawn. 

2. Derivation of matrix product states and realization of SU{2) invariance 

We are going to study quantum spin systems with local interactions. It is well-known 
that quantum system in d spatial dimensions can be mapped to classical systems 
in d + 1 dimensions. In this way, the ground-state properties of a quantum chain 
in the thermodynamical limit are encoded by a classical system on an unrestricted 
2-dimensional square lattice. Often, for numerical purposes, the quantum chain is 
mapped onto a kind of Ising model on a square lattice with chequerboard structure, 
for analytical purposes the mapping of quantum chains onto vertex models on periodic 
square lattices is more convenient. (The associated classical vertex model has nearest- 
neighbour couplings even for quantum spin-chains with interactions ranging farther than 
nearest-neighbours.) Here, the reasoning is based on equilibrium states, but obviously 
the derivation is more general and covers all steady state systems with local interactions. 
After mapping the quantum chain onto a vertex model, the correlation functions 
of the classical model on the full plane yield the correlations of the quantum chain, as 
is well-known. It is less well-known, but equally easy to understand that the partition 
function of a half-plane with arbitrary, but fixed boundary spins yields the coefficients of 
the ground-state of the quantum chain with respect to the standard basis. The evolution 
operator associated with a column of the full plane is known as the transfer matrix of 
the model. The corresponding objects of the half-plane are known as (lattice) vertex 
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Figure 1. a) Illustration of the relation of quantum spin chains and classical models 
on square lattices. The rows and columns of the square lattice define the action 
of the row-to-row and the column-to-column transfer matrix, b) The coordinates 
(...cri_iiTiCri+i...|'0) of the ground-state of the quantum system correspond to the 
partition functions of the classical model on the half-plane with fixed boundary spins 
...cTi-iiTiCri+i.... The column of the half-plane with fixed boundary defines the vertex 
operator. 



operators. With view to Fig. [T] these objects carry spin variables on the left, on the 
right and one spin variable on the top. When considering these objects as matrices 
where the spin variables on the left play the role of the row index, and the spin variables 
on right play the role of the column index, vertex operators are matrices with (local) 
spin state valued entries. The goal of the DMRG procedure may be understood as the 
computation of the optimal truncation of the infinite dimensional vertex operator to a 
finite dimensional matrix space. In this section we present the algebraic background for 
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the construction of general su{2) invariant MPS. 

We first summarize the basic algebraic constructions needed for our investigation of 
many-body quantum spin systems. We place particular emphasis on the compatibility 
with symmetry groups notably Lie groups. Eventually we will be interested in the Lie 
group SU{2) which is the reason for being specific from the beginning. 

We consider the class of matrix-product states 

|^) = Tr((7i-^2-...-^L) (1) 

where Qi is a square matrix with some auxiliary (index) space V and entries from a local 
quantum space V^ which we take as the ith copy of a su{2) spin-s space C^^"*"^. 

S'f/(2)-invariance of {ip) is guaranteed if for the representation Us of SU{2) in C^*"^^ 
there is a representation Uy in V such that Us applied to any element of the matrix g^, 
denoted by Usgi, yields a matrix identical to Uy^ ■ Qi ■ Uy where dots refer to matrix 
multiplication. Obviously we obtain with ([T]) 

Us(^Us® ... ® Usltp) = TiiUsgi ■ Usg2 ■ ■■■ ■ Usgi) 

= Tr(f/^i -g^-Uv-Uy^-g^-Uv-...- Uy' ■ gr. ■ Uv) 
= 1^) (2) 

The local condition for S'f/(2)-invariance can be written as 

Uv®Usg Uy^ = g, (3) 

meaning that the object g may be regarded as a tensor of the space V (g) C^''^^ (g) V* 
where V* is the dual space to V. (Note that the product A ■ B oi two linear maps A 
and B of the space V corresponds to the tensor product followed by a contraction of A 
and B viewed as tensors in V ®V*.) 

For the purpose of imposing discrete lattice symmetries like parity, i.e. invariance 
with respect to reflections, we adopt a different point of view. Let us consider tensors 
G from V ® C^^"*"^ ^V. As a local condition for S'f/(2)-invariance we demand 

Uv^Us®UvG = G, (4) 

and as a local condition for parity invariance we demand - as a sufficient condition - 
that G be symmetric with respect to exchange of "the first and the third index" when 
written in a canonical basis. 

The relation between g and G is realized by a 5'f/(2)-invariant tensor S from V ^V 
(and by the invariant dual tensor S* from V* Cg V"*). Concrete candidates for S (and 
S*) will be given shortly. 

The tensor 5* can equivalently be understood as a linear map from V* to V as the 
multiplication of an arbitrary element v of V* with S yields an object in V* ^ V ®V, 
the subsequent contraction over the first and second space yields an element v of V. 
Denoting S{v) := v we establish S as a map V* — )■ V. The S'f/(2)-invariance of S* as a 
tensor in V" V" is written as Uy Cg) Uy S = S from which we find 

S{vUy') = UyV = UyS{v). (5) 
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The object G is obtained by action of S on the third space of g 

G = id(»id®S g (6) 

which takes ([3]) into (j4j) thanks to ([5]). 

Conversely, we estabhsh S* as a hnear map V ^ V* with S'f/(2)-invariance 
S*{Uvv) = S*{v)Uy^ and for the concrete reahzations of S and S* we find S-S* = (-1)^-' 
and S* ■ S = (—1)^-' in spin-j subspaces. Hence, S and S* are invertible and 
g = id(g)id(g) S'^ G. 

Finally, we have to give explicit constructions for S (and S*) as S't/(2)-invariant 
states in V^® V". This is the only place where we make explicit use of the fact that our Lie 
symmetry group is SU{2). We take the space \^ as a direct sum of some irreducible spin- 
j representations where j = 0, |, 1, |, 2, ... . Each j may appear an arbitrary number of 
times, in which case we label the different orthogonal mupltiplets by an integer i. The 
space V is spanned by orthogonal states |(j, i),m) where the magnetic quantum number 
m varies from —j to +j in integer steps. SU{2) singlet states in V^ V and V* (g) V* 
are given by 

j,i rn=-j 
3,i ■rn=-j 

where 5"* is related to S by replacing the ket-states by the dual bra-states. 

Applying our above formulated definitions we find for v = {{j,i),m\ that S{v) = 
{—iy~"^\{j,i), —m). Conversely, for v = \{j,i),m) we have S*{v) = (— l)-'~™((j, i), —m\. 
Hence, the successive action of S and S* yields (—1)^"' id. 

3. Basic representation theoretical settings: 3j and 6j symbols 

Having spelled out the fundamental objects appearing as factors in S'f/(2)-invariant 
matrix-product states, the concrete calculations are straightforward. We want to use 
SU{2) singlets G in V^C'^^'^^^V. Having already allowed for reducible representations 
in V, we like to stress that for our applications we must deal with V as direct sum of 
more than one irreducible representation. This is so as for the most interesting case 
of s = I no singlet G exists if V is identical to just one spin-j multiplet. There is no 
half- integer spin in the tensor product decomposition of two spin-j 's! 

Let us consider in V (g) C^'^^"'^ ® V any spin multiplet {ji,ii) from the first factor 
space, the (only) spin multiplet J2 (= s) of the second space, and again any spin multiplet 
(J3, is) from the third factor space. Disregarding scalar factors, there is at most one way 
of coupling these multiplets to a singlet state. The coupling coefficients are known as 
3j symbols and the desired singlet is 

|(jl,«l),J2, (J3,^3)) : = 
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Figure 2. Introduction of ttie graphical notation for singlets appearing in product 
spaces. Three leg (vertex) symbols for |(ji, Ji), J2, Oa, «3)) and (Oi, Ji), J2, O3, is)! 
appearing in products of three irreducible representations. Two leg (edge) symbols 
for the singlets S resp. S** in T^ x T^ resp. V* xV*, and arcs for singlets in F x V^* for 
the same spin j . 
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Figure 3. Graphical notation for key arithmetic objects like Kronecker symbol and 
3j symbol. 
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r?ii,m2,»n3 



Jl J2 J3 

nil ^^2 ^^3 



\{ji,ii),mi) (g) \J2,m2) |(j3,i3),m3). 



Further below we will be using a graphical language for constructions and actual 
calculations. For instance, the singlet (|8]) and its dual are depicted by three straight 
lines carrying arrows, see Fig. [2J The coupling coefficients of the singlets shown in Fig. [2] 
are given in Fig. |3l More precisely, the objects shown Fig. |2]are obtained by multiplying 
the objects in Fig. [3] by states \{j,i),m) (or the dual) and summing over all magnetic 
quantum numbers m. 
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The singlet G can be written as superposition of these elementary singlets 

G= Y. ^:!!\Ul^^^^J2,{J3,^3)) (9) 

with suitable coefficients Al^'J^. Note that J2 has been suppressed as index- like argument 
of A as J2 is always identical to | and unique (for this reason no Z2 has been introduced 
above). Due to the symmetry of 3j symbols with respect to exchange of two columns 



J3 J2 Jl \ _ /_-.NJ^+j,-2+J3 





fiyi+n+]3 -J^ -J^ -J-^ fiQ\ 

m^ — — • I ™ ™ ™ / 

we conclude that 

^ii = i-^'^^^'^^'Ki:- (11) 

is a sufficient condition for parity invariance. Note that ji + J2 + js is always integer. 

Also note that only few combinations ji,J2,J3 need to be considered: if the triangle 
condition \ji — J2I ^ js ^ \ji + J2I or any condition obtained by permutations of the 
indices is violated, the three multiplets can not couple to a singlet. Since J2 = I, we are 
left with the combinations {jijjs} = {0, ^j, {|, 1}, {1, |}, .... This is a natural result, 
since ji and J2 = | only couple to ji ± ^. Let us denote by Uj the number of spin-j 
multiplets. By use of the symmetry (ITTl) we may reduce all possible coefficients A to a 
set of rij X rij+i/a matrices A-^'-^+Vz with matrix elements 

(a^'^+vA ,= A^^y^_ (12) 

3. 1 . Norm and transfer matrix 

Next, we want to calculate the norm {ip\il)) and the expectation value of the Hamiltonian 
{'4)\H\ip) in the thermodynamic limit. The computation leads to 

(^1^) = Tr(^+^i ■ g^g2 ■ ... ■ gtdL), (13) 

where 5-+ G V* O (C2"+^)* ®V is the dual oi g ^V ® C^'"^^ ® V* and the contraction 
over the second space is implicitly understood in g^g. Hence T := g^g is a linear map 
V ®V* — 7- V ®V* . For the computation of the norm we employ the transfer matrix 
trick yielding for the r.h.s. of ( IT3|) 

(^|^) = Tr(T_^^) = 5^A^ (14) 

L times ^ 

where the sum is over all eigenvalues A of T. Obviously, in the thermodynamic limit 
only the largest eigenvalue (s) contribute. 

The computation of the leading eigenvalue is facilitated by the singlet nature of the 
leading eigenstate. There are not many independent singlet states vnV ®V* . A (j, -j) 
multiplet in V and a (j, i) multiplet in V* couple to a singlet iff j = j (with arbitrary i 
and i). The (normalized) singlet is given by 

1 ^ 

y/2j + 1 ^ . 
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Figure 4. Illustration of identity (|18p : The summation over two internal magnetic 
quantum numbers leads to an "arc" -singlet. 



Graphicially, this singlet is depicted by a link carrying an arrow pointing from the V* 
to the V space. 

The action of the transfer matrix T onto a singlet a{j] i, i) produces similar singlets 
where j is changed by ±i. The a{ji] ii,ii) — cr{jr; ir,ir) matrix element of T is 

where ji = jr ± |. (Expression flT6|) is the analogue of eq. (6) in [2], but here we do 
not impose condition (10) of ^.) The matrix T has a simple block structure with zero 
diagonal blocks and non-zero secondary diagonal blocks. From (TTT!) we conclude that 
the matrix is symmetric. The defining blocks are 

In expression (fT6ll the coefficients A and the complex conjugate A* derive from the 
explicit appearance in (|9]) and the prefactor is obtained from the identity of 3j symbols 
(if ji is found in the product of J2 and J3) 



y '] '" 'M '' '" ■" =^— 5,./5™,„.', (18) 

^-^ \ m[ m2 m^ j \ nii 7/12 m^ j 2ji + 1 ' ^ ' ^ 

illustrated in Fig. HJ The computation of the matrix elements of T is graphically 
presented in Fig. |5l 

For a given space V with a certain number of spin-0, |, 1, | ... multiplets and a 
certain set of parameters A{ji,ii]jr,ir), the transfer matrix has to be diagonalized in 
the singlet space. The total number of non-zero coefficients A{ji, ii]jr, ir) is J2j njnj+1/2 
(with sum over j = 0,|,1,|,... and Uj denoting the number of spin-j multiplets). The 
total dimension oi V ^V* is { ^ • nj{2j + 1)) , but the singlet subspace is much lower 
dimensional: X^i^?- -Due to the still high dimensionality, the diagonalisation in the 
singlet space has to be done numerically. The eigenvalues come in pairs ±A. The 
two largest eigenvalues and the corresponding eigenstates determine the physics in the 
thermodynamic limit. 
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Figure 5. Depiction of the transfer matrix T in terms of local vertices and the action 
onto singlet states. 



3.2. Nearest-neighbour couplings 

We are interested in the spin-| Heisenberg chain with nearest-neighbour interaction 
with Hamiltonian 

L 

H = Y,SiSi+i. (19) 

1=1 

The local Hamiltonian is /iz,z+i = SiSi^i = 1/4 — P^" where Pq" is the projector 
onto the nearest-neighbour singlet space. We want to determine the matrix-product 
state with minimal expectation value of the total Hamiltonian H = Yli^i- -Due to 
translational invariance this is achieved by minimizing the expectation value of a single 
local interaction. In analogy to ( 1T^ we obtain 

{iP\Pr\i^) = Tr(T T__^) = Ao^-2(0|T2|0), (20) 

L-2 times 
where we assumed P^"" to act on sites 1 and 2. T2 is a modified transfer matrix acting 
in V ^ V*, |0) is the (normalized) leading eigenstate of the transfer matrix T and we 
kept the only term dominating in the thermodynamical limit. Hence 



|pnn| 



ii^m 



Ao'{0\T,\0). 



(21) 



The computation of the matrix elements of T2 is described graphically in Figs. |6] and 
[71 In contrast to the transfer matrix T the modified matrix T2 is block diagonal with 

a{j;ii,ii) — cr(j; v, ir) matrix element 

{j;tiMT,\j;t^X) = 77T-^5^lfei)*- (22) 



(2j + l)2 ^''^A n^^r 



The matrices B are given by 



(23) 



ji 
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Figure 6. Depiction of tlic modified matrix T2 in terms of local vertices. 
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Figure 7. The summation over two internal magnetic quantum numbers leads to an 
"edge" -singlet. 
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Figure 8. Depiction of the modified matrix T3 in terms of local vertices. 



where only the values ji = j ± | lead to non-zero terms. (Equations fl22|23l) are the 
analogue of eq. (26) in 12].) Using this and the (sufficient) condition flTTl) for parity 
invariance we find 



B 



<3,3 



^i-V2.i y^i-V2.i _|_ ^ij+V2 ^i.i+Va 



(24) 



3.3. Next-Nearest-neighbour couplings 



The next-nearest neighbour interactions are manageable, too. In the thermodynamical 
limit we find 



|pnnn|, 



Ao-^(0|T3|0). 



(25) 
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Figure 9. The summation over three internal magnetic quantum numbers leads to a 
6j symbol times an elementary vertex. 



The computation of the matrix elements of T3 is described graphically in Fig. [8] and 
Fig. ini In contrast to the modified transfer matrix T2, but like T, the matrix T3 has zero 
diagonal blocks and non-zero secondary diagonal blocks. The a{ji]ii,ii) — a{jr]ir,ir) 
matrix element of T3 is 



{jl]il,k\T3\jr;ir,ir) 



1 



v/(2j, + l)(2j;, + l 






where ji = jr ± \ and C^''^'' is given by 



(JjlJr _ y^C_]^yV2+Jl+i3 



J1J3 



jl V2 Jr 
J3 V2 Jl 



j\Jl,nj\Jlj3j\j3,Jr 



(26) 



(27) 



(There is no analogue to eqs. 026II27P in [Il[2]-) Only three combinations of ji,J3 yield 
non- vanishing contributions. For jV = ji + 1/2 only (ji, J3) = {ji - 1/2, jz), {ji + 1/2, ji), 
{jl + 1/2, jl + 1) are relevant. From this and the (sufficient) condition (ITT]) for parity 
invariance we find 

j 1/2 J + V2 W^j-V2.A'^^i-V2J^iJ+V2 



QJd+^/2 



+ 



3 V2 3 - V2 



J V2 J + V2 
J V2 J + V- 



j + 1 1/2 J + V^ 
where the 6j symbols evaluate to 

j 1/2 J + V2 1 (-1)^^+^ 
J V2 3 - V2 



■^' V2 J + 72 I ^iJ + l/2^j + l/2j + l /^i + l/2,i + iy 



(28) 



J V2 3 + V2 

J V2 3 + V2 

j 1/2 3 + V2 

j + 1 1/2 J + V2 



2j + l 

(-l)^-^+i 
(2j + l)(2j + 2)' 

(-1)^^' 
2j + 2' 



(29) 






-0.283412 














0.183208 














0.679798 



Diagrammatics for SU{2) invariant matrix product states 13 

with an exception for j = where the first hsted 6j symbol has to be taken as 0. 

4. Results 

For the nearest-neighbour spin-^ Heisenberg chain we found that already a few low- 
dimensional multiplets in the matrix space V yield excellent results, e.g. the ground- 
state energy differs from the exact result cq = 1/4 — In 2 by about 5 ■ 10~^ and dimer 
correlations are off the exact results by about 10~^. This is achieved with uq = 4, 
ni/2 = 4, rii = 3, ria/a = 2, n2 = 1 (^^=0 for s > 2). Note that the attempt to include 
higher spin multiplets at the expense of reducing the low-spin multiplets is not successful 
as for instance Ug = 1 for all s leads to the simple dimer (Majumdar-Ghosh) state. 

The minimization of the energy expectation value yields the following list of 
coefficients 

" 1.0 



^0,1/2 



0.737765 

0.257549 -0.860519 

0.610288 -0.093768 0.050733 

0.191246 -0.631258 0.181757 

0.510676 

-0.5075510 0.714082 

-0.292858 0.712182 

and a value of the ground-state energy which compares well with the exact value ^ 

gnum ^ _o.443 092 175..., e^,^^^* = -0.443 147 180... (31) 

Note that A^'^^^ is strictly diagonal, and the other matrices have strictly zero entries 
above the diagonal due to a "gauge" freedom. The MPS is invariant under a 
transformation A^''^'' — )■ 0~^A^''^''Oj^ where Oj are arbitrary orthogonal n^ xn^-matrices. 
Also note that the entries of the matrices are strictly real. 

Next we give numerical results for the dimer-dimer correlation function 

Dn:={{S,-S2){Sn+l-Sn+2)). (32) 

For n — 7- oo this converges to Cg, hence it is more instructive to study the connected 
dimer correlations 



^VM 



^1,3/2 



^3/.,2 



0.265736 
-0.607939 



(30) 



Dr 



-^n J-^OQ- 



(33) 



For n = 


= 2,3, 


4 the values are known exactly [21 


MPS calculations 










^num 


= +0.060 639..., 


^exact 






r-)num 
^3 


= -0.027 838..., 


r-) exact 
^3 






jjnnm 


= +0.018 986..., 


r-j exact 
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which we use for comparison with our 

= +0.060 824..., 

= -0.027 737..., (34) 

= +0.018 928.... 

We expect that the absolute numerical accuracy is similar also for larger distances n 
of the local dimer operators. These results are plotted in Fig. [10] as (— l)"!)^ versus 
n. Note that all (— 1)"'D„ are positive which implies a sublattice structure with sign 
alternation of the correlations. There is no long-range order for the spin-| Heisenberg 
chain with nearest-neighbour interactions. 

Next we are interested in the frustrated spin-| Heisenberg chain with Hamiltonian 

L 

H = J2 {SiSi+i + aSiSi+2) . (35) 

1=1 
The system shows algebraically decaying dimer-dimer correlations for a < 0.2411..., see 
Fig. [HI and long-range dimer order for a > 0.2411..., see Fig. [TO]b). The critical value 
ttc = 0.2411 was established in [361 EZ]- The dimer-dimer correlations are fitted well by 
algebraic curves for all a < 0.2411.. with a-dependent exponent. From field theoretical 
considerations the exponent is expected to be identical to 1. We attribute the deviations 
to logarithmic corrections for a < 0.2411.. which apparently vanish at ac- 

As an illustration of spin-spin-correlation functions we present results for the 
nearest- and next-nearest-neighbour cases in table [H For a = the exact values are 
known [3l] Cf'^"* = eg^'^"*, see (E]), and [35] Cf^^^* = 0.182 039.... Even for C3 our 
numerical value deviates from the exact value less than 10~^. Due to the variational 
nature of our calculations the expectation values of the energy are strict upper bounds 
for the ground-state energy. Note that for cases a = 0.1, 0.2411, 0.3, 0.4 our results are 
lower than those given in [32] with small differences ranging from 10~^ to 10"'^. Hence, 
the deviation of these DMRG results from the true ground-state energy must be of the 
same order or even larger. In DMRG calculations there are two sources of errors: (i) 
truncation of the Hilbert space and (ii) finite-size effects due to the finite length of the 
considered chains. In our approach we deal with the strictly infinitely long chain. The 
numerical computations of the data presented in this section were done by use of Maple 
13 on a laptop computer. The total computation for the seven cases of the frustration 
parameter took about 1 hour. 

A more complete study of the correlation functions and of the physics of frustrated 
systems with a > 0.5 will be presented elsewhere. Here we like to note that for a > 0.5 
the matrices replacing ( l30ll will contain intrinsically complex numbers. 

5. Conclusion 

We showed how to employ systematically su{2) invariance for matrix product states and 
how to carry out the variational computation of the ground state energy in a numerically 
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Figure 10. Dimer-dimer correlations of the spin-i Heisenberg chain with frustration 
parameter a: Linear plots of (— 1)"Z?„ for a) a = 0,0.1,0.2,0.2411 and h) a = 
0.2411,0.3,0.4,0.5. Note the long-range order in plot b) for a > 0.2411. 



a 



0.0 



0.1 



0.2 



etc. 



0.3 



0.4 



0.5 



C2 


-0.443092 


-0.442655 


-0.440916 


-0.439574 


-0.436475 


-0.420659 


-0.375000 


C-s 


0.181942 


0.173570 


0.162233 


0.156176 


0.144794 


0.100870 


0.0 


eo 


-0.443092 


-0.425298 


-0.408469 


-0.401920 


-0.393037 


-0.380311 


-0.375000 



Table 1. Numerical values for the correlations C„ := (SiSn+i) for n = 1,2 and 
various values of the frustration parameter a. The ground-state energy cq = Ci + aC2 
is given in the last row. 
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Figure 11. Dimer-dimer correlations of the spin-i Heisenberg chain with frustration 
parameter a < 0.2411: Double log plot of (— 1)"Z3„ versus distance n and algebraic 
lines c/n^ from fits to the data for n in the range [10, 25]. 



most efficient manner. The algebraic computations for the sM(2)-invariant bulding 
blocks were put in diagrammatic formulation. As an example we used the (frustrated) 
spin- 1 Heisenberg chain with nearest- and next-nearest neighbour interaction. Our 
algebraic constructions led to the main results f lTB]) for the transfer matrix, and (12211231) 
and (I26II27P for the modified transfer matrices, where the coefficient matrix has to satisfy 
the relation (llip for parity invariance. 

Our calculations are very similar to those of [H |2] who applied the method to 
gapped spin-1 chains and spin ladders. The variational MPS calculations for the spin-| 
Heisenberg chain are demanding on their own: the model shows algebraically decaying 
correlation functions and it is by no means clear if signatures of this decay can already 
be seen in variational MPS calculations with only few multiplets in the matrix space. 
Also, in contrast to [1] [2], the matrix space we had to deal with consists of all integer 
and half-odd integer spin multiplets. In our concrete calculation we used a matrix 
space composed of 4 singlets, 4 doublets, 3 triplets, 2 quadruplets, and 1 quintuplets. 
We managed to calculate the ground-state energy within a precision better than 10~^. 
Also, the correlation functions were computed within an accuracy of the order 10~^ and 
allowed for the identification of the scaling dimension in the case of critical frustration 
(a = 0.2411). 

The actual numerical calculations like matrix diagonalizations were reduced by the 
algebraic su{2) implementation from a 1156-dimensional to a 46-dimensional space. The 
systematic inclusion of more and higher spin multiplets is obvious. Generalizations of 
these calculations are straight forward, e.g. to spin-S* Heisenberg chains with competing 
interactions. We are convinced that a systematic application of symmetries to the MPS 
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analysis of quantum spin chains will provide high quality data with only small truncation 
errors. 
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